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Hydriding properties of uranium alloys have been studied to search for new hydrogen storage materials to be 
applied to hydrogen energy systems. Application of uranium-base hydrogen storage materials can be expected 
to alleviate the risk, as well as to reduce the cost incurred by globally-stored large amounts of depleted uranium 
left after uranium enrichment. Various uranium alloys have been examined in terms of hydrogen absorption- 
desorption properties, among which UNiAI intermetallic compound showed promising characteristics, such 
as lower absorption-desorption temperatures and better anti-powdering strength. First principle calculation 
has been carried out on UNiAI hydride to predict the change of crystal structure and the lattice constant with 
increasing hydrogen content, which showed this calculation to be promising in predicting candidates for good 


hydrogen absorbers. 


Keywords: Depleted uranium, Hydrogen storage material, Uranium alloys, UNiAI intermetallic compound, First principles 


calculation 


DOI: 10.13538/j.1001-8042/nst.26.S 10312 


I. INTRODUCTION 


The hydrogen energy system requires a large-scale station- 
ary storage facility at hydrogen generation plants. Common 
forms that hydrogen can be stored in include compressed gas, 
cryogenic liquid, and metal hydrides [1]. Each of these stor- 
age methods has advantages and disadvantages, as well as re- 
quires some degree of energy consumption, that mainly aris- 
es in gas compression, maintaining cryogenic temperature, or 
heating metal hydrides for desorption, respectively. One of 
the advantages of using metal hydrides is the inherent safety. 
It only requires moderate pressure when storing and passive- 
ly contains hydrogen in case of power loss. One disadvantage 
of metal hydrides for mobile storage is limited mass energy 
density. This is not a practical concern to the stationary sys- 
tem, or rather, they may produce as high a volume energy 
density as other methods. The cost of the metal hydride s- 
torage system would largely depend on the material and the 
operational temperature range. Yamawaki et al. have previ- 
ously proposed a surplus electricity storage system utilizing 
depleted uranium as a base storage material [2, 3]. 

Depleted uranium has been produced abundantly as waste 
from the uranium enrichment process but has never been used 
effectively. It is usually stored idly as uranium fluoride gas or 
in solid form. On the contrary, uranium is one of the high- 
est capacity hydrogen storage materials. It readily absorbs 
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hydrogen to form uranium hydride, UH;, at room tempera- 
ture and discharges hydrogen reversibly at an elevated tem- 
perature. In fact, urantum has commonly been used to store 
tritium, the radioactive isotope of hydrogen, since the room 
temperature desorption pressure is low enough to securely 
contain it. On the other hand, discharging hydrogen at a dis- 
sociation pressure of 10° Pa requires heating the hydride up to 
about 700K. The discharging temperature is rather high for 
hydrogen storage applications incurring significantly large 
operational energy loss. Further, uranium undergoes large 
volumetric expansion upon hydrogenation due to the large 
change in crystal structure, which causes severe powdering of 
the material leading to low thermal conductivity and increas- 
ing pyrophilicity. Hydrogen absorption-desorption properties 
of various uranium alloys have been studied in order to find 
alloys with improved anti-powdering capacity as well as high 
desorption pressure by Yamawaki and co-workers [3—10]. 

The objective of this paper is to summarize the status and 
approach development of uranium alloys as hydrogen storage 
material. 


Il. HYDROGEN ABSORPTION-DESORPTION 
PROPERTIES OF URANIUM ALLOYS 


A. U-Zr and U-Ti alloys 


Hydrogen absorption-desorption properties and the reactiv- 
ity with air of U-Zr and U-Ti alloys were investigated for use 
in depleted uranium based energy storage systems as well as 
for tritium storage systems [3-8]. Some of the U-Zr and U- 
Ti alloys have hydrogen absorption-desorption properties as 
good as those of uranium. Further, the pyrophilicity of Zr- or 
Ti-rich uranium alloys is more moderate than that of U-rich 
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TABLE 1. Hydrogen capacity and powdering behavior of U-Zr and 
U-Ti alloys [3] 


H capacity [H]/[UM.] Powdering after cycling 


Alloy Maximum Rechargeable Diameter Abs.-Des. 
(um) Cycles 

U (ref.) 3 3 20-30 2-3 
UZrn 4 3 = 19.5 
UZr 5 2.7 100 11 
UZr, 6.5 1-3 1—400 12.5 
UTi,p 3.7 2.7-3.7 =e 13 
UTi 5 2.3 —4.5° -* 14.5 
UTi, 7 1.2 - 5.5° 1-10 12 


* The higher capacities with utilizing low UTi H+ phase plateaus; 
P Fine powder except a few granules. 


alloys due to the higher resistance against powdering upon 
hydrogenation [3, 7]. Table 1 summarizes the hydrogen stor- 
age capacities of U-Zr and U-Ti alloys [3]. On hydrogenation, 
U-Zr and U-Ti alloys decompose into UH; and ZrH, (x < 2) 
and UH; and UTi Hy (x < 5) phases, respectively, so that 
P-c-T diagrams consist of the plateau regions of correspond- 
ing phases. Because of low equilibrium H pressure of ZrH,, 
phases, the rechargeable hydrogen capacities of U-Zr alloys 
depended on the fraction of UH, phase (or U content in the 
alloy), except UZr, where the UH, region becomes smaller 
after cycles. In U-Ti alloys the ternary phase UTi,H, for- 
mation consumes some uranium to narrow the UH; plateau 
region, while utilizing the lower pressure UTi,H; plateau can 
extend the rechargeable capacity. Alloying U with Zr or Ti 
basically keeps the equilibrium pressure of UH3, while pow- 
dering is greatly reduced, especially with high content of Zr 
or Ti. Table 1 also shows the degree of powdering after up to 
=~ 20 hydrogen absorption-desorption cycles [3]. UZr, UZr,, 
and UTi, alloys showed significant resistance against severe 
disintegration into fine powders even after more than 10 cy- 
cles. UZr especially possess a good balance of large capacity 
and high resistance against powdering. 


B. U-Mn, U-Fe and U-Ni alloys 


Uranium-rich intermetallic compounds, namely UM (M 
= Fe, Mn, Ni), commonly forms in these binary systems. 
Cubic MgCu,-type UM, phase is the only other intermetallic 
phase in the U-Mn and U-Fe systems, while four additional 
phases form in the U-Ni system, including UNi; and UNi, 
of known AuBe,;- and Cu,Mg-type crystal structure, respec- 
tively, as well as U;Ni, and U7Nig. The hydrogen absorp- 
tion properties of these compounds have been studied [6, 10], 
where none of the alloys except the U;M-type compounds 
were found to be reactive to hydrogen. All UM type com- 
pounds absorbed hydrogen to form UH, and the correspond- 
ing UM, phases. The P-c-T isotherms of these alloys basical- 
ly consist of a wide single plateau of UH. As shown in Fig. 1, 
the plateau pressures in these U,M alloys are all higher than 
pure U. Among the alloys, U,Fe has the highest desorption 
pressure. 
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Fig. 1. (Color online) Hydrogen desorption pressure of UsM type 
alloys in comparison with uranium. 


C. UCo 


UCo, one of the six known intermetallic compounds in 
the U-Co system, is known to form UCoH,, [9, 11]. Our 
study on the hydrogen absorption—desorption behaviors of 
the compound has shown that there are two stages of desorp- 
tion reactions at higher temperatures than 473 K, while only 
a single stage was observed below 423K [9]. UCo readi- 
ly absorbs hydrogen to form UCoH, 7 at temperature below 
423 K under the hydrogen pressure of 10° Pa. It simply dehy- 
drogenates reversibly at 423 K, while two-phases of UH, and 
UCo,, appear during the desorption of UCoH, 7 above 473 K. 
This is due to the following two competing reactions: 


UCoH, 7 = 0.5UCo, + 0.5UH;+0.6H,, (1) 


UCoH, 7 = UCo + 1.35H3. (2) 


The standard Gibbs energy changes, AG? and AG$, for 
reactions (1) and (2), respectively, as functions of T and hy- 
drogen pressure, p, were estimated as follows 


AG® = 8.4 x 10*—174T + 0.6RT In(p/p,), (3) 


AG? = 7.4 x 104 — 147T +1.35RT In(p/p,), (4) 


here, po is the standard pressure (101325 Pa), R is the gas 
constant (8.3 J/(mol K)), and each Gibbs energy is calculated 
per mole of UCoH, ;. Fig. 2 plots the AG? and AG$ vs. T 
in lines with unfilled and filled symbols, respectively, for two 
cases of p. Solid lines with circles are for p = 104 Pa, while 
broken lines with squares are for p = 10° Pa. For lower p 
(= 10* Pa), AGS is lower than AG? over a wide temperature 
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Fig. 2. (Color online) (a) Standard Gibbs energy changes of two 
competing reactions on dehydrogenation of UCoH, for the hydro- 
gen pressure p of 10* Pa (Solid lines) and 10° Pa (dotted lines). The 
reaction (1) to decompose into two phases and (2) to reverse to UCo 
are shown in the lines with unfilled and filled symbols, respective- 
ly; (b) Equilibrium desorption pressures based on the enthalpy and 
entropy changes given in Ref. [9]. 


range below 650K, while at p = 10° Pa, AG9 is higher at 
T > 400K. This results in the observation of a single stage 
reaction (2) at lower T and lower p, while two stages, includ- 
ing the decomposition into UCo, and UH, result at higher T 
and higher p. 

While desorption operations may be limited in the p range, 
UCo can reversibly absorb and desorb hydrogen at a rela- 
tively low temperature of 423 K with the desorption pressure 
significantly higher than that of UH}. 


Nucl. Sci. Tech. 26, S10312 (2015) 


D. UNiAI and other Fe,P-type uranium compounds 


UNiAI intermetallic compounds have a hexagonal Fe,P- 
type (or also known as ZrNiAl-type) crystal structure. As 
illustrated in Fig. 3(a), the unit cell consists of a stack of t- 
wo kinds of planes, one with three U and one Ni atom and 
the other with three Al and two Ni atoms. Table 2 shows 
more details of the atomic positions in terms of Wyckoff’s 
notations and representative partial coordinates. In the origi- 
nal crystal structure, there are two U atom triangles per unit 
cell, which with two Ni atoms on the adjacent planes make 
tetrahedral interstitial sites sharing the triangle, as illustrat- 
ed also in Fig. 3(a). At hydrogen pressures below 10° Pa, 
UNiAI absorbs up to two hydrogen atoms per unit cell, or 
around 0.7 H per UNiAI formula unit (F.U.) (H/F.U.) [11- 
13]. The neutron diffraction study of a deuteride, UNI AID} 7, 
showed that D atoms occupy the U3Ni,-type interstitial sites. 
At a higher hydrogen pressure it can absorb up to around 2 
H/EU. [14]. Yamamoto et al. showed that the crystal struc- 
ture of UNiAIH,, or UNiAID,, has significantly different 
atomic positions while retaining basic hexagonal structure. 
Table 2 also shows the detailed crystal structure of the deu- 
terides, including the atomic positions. In UNiAID, , forma- 
tion the uranium positions shift to increase x in the repre- 
sentative coordinate, (x, 0, 1/2), of 3g position. Aluminum 
atoms also move within the plane in the same direction as 
that of U atoms, while Ni atoms in the 1b position moves a- 
long the c-axis to z = 0 from z = 1/2. These changes are 
illustrated in Fig. 3(b), 3(c), and 3(d) show U and Ni atom 
positions on the z = 1/2 plane for UNiAl, UNiAIDo,, and 
UNiAID, 5, respectively. U atoms shift to increase x from 
0.580 to 0.592 and 0.659 makes the U triangle, forming an 
U;Ni,-type interstitial site larger by 4% and 17%, respective- 
ly, which presumably is required to secure the U-D distance. 
It is noticeable that these shifts also make the other U triangles 
similar equilateral ones. Deuterium occupying sites are also 
shown in Fig. 3(c) and 3(d) with three different triangle sym- 
bols for U;Nij-type D(1), U;Al,-type D(2), and U,Nij-type 
D(3) sites, respectively. The final D site structure suggest- 
s that it is the atomic position shifts that really enabled this 
compound to accept around 6 hydrogen atoms per unit cell. 

Hydrogen desorption isotherms of UNiAI have been stud- 
ied in detail by Yamanaka et al. [14]. Fig. 4 shows the hy- 
drogen desorption pressures at two plateau regions, between 
«x and B and B and y phases, corresponding to the following 
nominal reactions, respectively: 


UNiA1 + 1/3H, = UNiATH, 66; (5) 


UNiAIHg 65 + 2/3H, = UNiAlHy, (6) 


here, H compositions are all nominal values based on the 
number of H atoms occupying the 4h site and the 4h, 3g, and 
2e sites for B and y phases, respectively. The figure indi- 
cates that the desorption pressure for the reaction (6) reaches 
around 10° Pa at 400 K, which is sufficiently high for hydro- 
gen storage use of the material. The lower plateau reaches 
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(b) 


Fig. 3. (Color online) (a) Crystal structure of UNiAl showing four U,Ni'-type tetrahedral sites; (b) atomic arrangements on the U plane at 


z = 1/2 in (a) UNiAI, (c) UNiAIDo7 and (d) UNiAID, 9. 


10° Pa at around 550K, which may be acceptable, while us- 
ing only higher plateau region allows recycling of around 2/3 
of the total hydrogen capacity. 


II. FIRST PRINCIPLES CALCULATIONS BASED 
APPROACH TO NEW URANIUM HYDRIDES 


There are around 20 known ternary uranium compounds 
of ZrNiAl-type crystal structure. A few of them, as well 
as their quaternary mixtures have been investigated in terms 
of hydrogen absorption properties [16-19]. Of all those that 
have been examined, UNiAI has the largest hydrogen capaci- 
ty at around 2H/F.U., while that all others hold no more than 
around 1H/F.U. This suggests that the relocations of metallic 
atoms described in the previous section are of great impor- 
tance to the large hydrogen capacity. In searching for other 
ZrNiAl-type U intermetallic compounds that can hold large 
amount of hydrogen, knowing the possibility of such atomic 
relocations would provide a very helpful guidance. First prin- 
ciple calculations of potential hydrides can be used to exam- 
ine that possibility. To start off we have carried out the cal- 
culations of UNiAlH,, with n up to 2, which shows that the 
crystal unit shown in Fig. 3(a) holds up to 6 H atoms, to learn 
how the theory can inform us of the potential hydrogenation. 
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Fig. 4. (Color online) Hydrogen desorption pressures of UNiAI re- 
produced from P-x-T diagram in Ref. [14]. 
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TABLE 2. Crystal structures of UNiA] and its deuterides 


UNiAI [11] UNiAID 7 [10] UNiAID,, [11] 
Phase in Ref. [14] a y 
DFT model Structure A Structure B 
in Section III 
Lattice 0.6732 0.6962 0.7171 
constants c 0.4036 0.3992 0.3976 
Wyckoff’s for Coordinates & occupancy 
sg 189, P62m 
U 3g (0.580, 0, 1/2) 1 (0.592, 0, 1/2) 1 (0.659, 0, 1/2) 1 
Ni (1) lb (0, 0, 1/2) 1 (0, 0, 1/2) 1 — 
Ni (1’) la — — (0, 0, 0) 1 
Ni (2) 2c (1/3, 2/3, 0) 1 (1/3, 2/3, 0) 1 (1/3, 2/3, 0) 1 
Al 3f (0.219, 0, 0) 1 (0.239, 0, 0) 1 (0.347, 0, 0) 1 
D (1) 4h — (1/3, 2/3, 0.488) 1/2 (1/3, 2/3, 0.436) 0.40 
D (2) 3g — — (0.331, 0, 1/2) 1 
D (3) 2e — — (0, 0, 0.43) 1/2 
D (4) 6j — — (0.207, 0.233, 0) 0.09 
D (5) 3f — — (0.577, 0, 0) 0.15 
-©-la (exp.) © le (exp.) 
—@- la (DFT @ Ic DFT 
-83 1 T T ] 0.72 = i ) 1 D i 
-F Structure A | UNiAIH 
—@- Structure B 0.71 L n 
S -83.5 L J 
o (a) 07 | 4 0.42 
M -84 DFT 
z 3 0.69 L Ej 
D & 40.418 
& -845l = 0.68 + 2 
= 
= 0.67 
= 85 [ ` f 4 0.4 
H -i H . 
3(UNIAIH „) F 3(1-n/2H, 0.66 L i 
[ a 
-85.5 | ! 0.65 0.39 
0 0.5 1 15 2 -0.5 0 0.5 1 1.5 2 2.5 
n =H/UNiAl n=H/UNiAI 


Fig. 5. (Color online) (a) Total enthalpy of 3(UNiAIH,,) + 3(1 — n /2)H, system from the DFT calculations as a function of n in two metallic 
atom structures; (b) the lattice constants, la and Ic, as a function of H concentration, n (= H/UNiAI ) for both of DFT and experiments. 


The first-principle calculations based on the density- 
functional theory (DFT) are performed using the Vienna 
ab-initio simulation package, VASP, with generalized gra- 
dient approximation and projected augment wave potentials 
[20, 21]. We employed a hexagonal unit-cell containing three 
UNiAI formula units as shown in Fig. 3(a) with its Brillouin 
zone with 3 3 5 k-point sampling using the Monkhorst-Pack 
scheme [22]. The plane-wave energy cutoff is 350eV. Atom- 
ic positions and cell size relaxations are iterated until the total 
energy between two ionic steps is smaller than 0.0005 eV. 

Two structures of metallic atoms were considered in the 
calculations. Structure A is based on the original UNiAl 
structure, while Structure B is based on the experimentally 
observed UNiAID, , crystal structure. Key difference is that 
1/3 of Ni atoms locate in 1(b) (0, 0, 1/2) positions in Struc- 


ture A but in 1(a) (0, 0, 0) in Structure B. Hydrogen atoms 
are initially placed in positions where D atoms were found to 
occupy in UNiAID,7 and UNiAID, 5. Details of calculation 
will be reported in the future, while the preliminary results of 
DFT calculations have supported various observations in the 
experiments. For example, Fig. 5(a) shows the total enthalpy 
of the 3(UNiAIH,,) + 3(1—n/2)H, system from the DFT cal- 
culations as a function of n in two metallic atom structures. 
Structure A has lower energy for n < 2/3 (2 H atoms in the 
crystal unit shown in Fig. 3(a)), while for n > 1, structure B 
becomes more energy favored. The most stable crystal struc- 
ture in DFT calculations also supports the observed shifts in 
atomic positions of U and Al along the z-axis (listed in Ta- 
ble 2). Fig. 5(b) plots the lattice constants, la and Ic, as a 
function of H concentration, n (= H/UNiAI) for both DFT 
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and experiments. The characteristic trends of the changes in 
the lattice constants show an increas in la while decreasing 
in lc with increasing n, are qualitatively in good agreement 
between the calculation and experiments. Although these are 
still preliminary results, as quantitative comparisons between 
calculations and experiments need to be examined in more 
detail, the results at least qualitatively support the experimen- 
tal trends suggesting that the approach would be able to give 
a list of candidates for good hydrogen absorbers. 


IV. CONCLUSION 


A hydrogen energy age is coming soon, where metal hy- 
dride will be used as a hydrogen storage material. The al- 
loy of uranium with Zr, Ti, Mn, Fe, Ni, and Co each has 
been examined in terms of its potential to be used as a hy- 
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The results have supported the experimental trends, suggest- 
ing the present approach to be promising in giving a list of 
candidates for good hydrogen absorbers. 


ACKNOWLEDGEMENTS 


The authors thank Mr. Kunihiro Ito of the Nuclear De- 
velopment Company for supplying uranium sample used to 
produce the UNiAI specimen of this study. 


[12] Yamamoto T, Kayano H, Yamawaki M. Crystal structure of hy- 

dride and deuteride of unial. J Alloy Compd, 1994, 213-214: 

533-535. DOI: 10.1016/0925-8388(94)9098 1-4 

[13] Yamamoto T, Ishii Y, Kayano H. Deuterium absorption proper- 

ties and crystal structure of UNiAl. J Alloy Compd, 1998, 269: 

162-165. DOI: 10.1016/S0925-8388(98)00221-7 

[14] Yamanaka S, Iguchi T, Fujita Y, et al. Study on the hydrogen 

solubility in UNiAl. J Alloy Compd, 1999, 293-295: 52-56. 

DOI: 10.1016/S0925-8388(99)00300-X 

[15] Drulis H, Petrynski W, Stalinski B, et al. Hydrogen absorp- 
tion properties of selected uranium intermetallic compounds. 
J Less-Common Met, 1982, 83: 87-93. DOI: 10.1016/0022- 
5088(82)90172-2 

[16] Jacob I, Hadari Z, Reilly J J. Hydrogen absorption in ANiAl (A 
= Zr, Y, U). J Less-Common Met, 1984, 103: 123-127. DOI: 
10.1016/0022-5088(84)90366-7 

[17] Raj P, Sathyamoorthy A, Shashikala K, et al. The hydriding 
behaviour of U(Fe1_~ Niz) Al system (0 < x < 0.75) and 12x 
x magnetic studies on U(Fe1_~ Nis )AlHo.s. J Alloy Compd, 
2000, 296: 20-26. DOI: 10.1016/S0925-8388(99)00525-3 

[18] Maskova S, Havela L, Danis S, et al. Onset of magnetism s- 

timulated by H absorption in UNiAl-UFeAl. J Phys Conf Ser, 

2011, 303: 012011. DOI: 10.1088/1742-6596/303/1/012011 

[19] Raj P, Sathyamoorthy A, Dhar S K, et al. Effect of deuterium 

absorption on the heavy-fermion compound UPdIn. Phys Rev 

B, 2002, 66, 214420. DOI: 10.1103/PhysRevB.66.214420 

[20] Kresse G and Hafner J. Ab initio molecular dynamics for liquid 

metals. Phys Rev B, 1993, 47: 558-561. DOI: 10.1103/Phys- 

RevB.66.214420 

[21] Kresse G and Furthmueller J. Efficient iterative schemes for ab 
initio total-energy calculations using a plane-wave basis set. 
Phys Rev B, 1996, 54: 11169-11186. DOI: 10.1103/Phys- 
RevB.54.11169 

[22] Monkhorst H J and Pack J D. Special points for Brillouin- 
zone integrations. Phys Rev B, 1976, 13: 5188-5192. DOI: 
10.1103/PhysRevB.13.5188 


$10312-6 


